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We propose qubits based on shallow donor electron spins in germanium structures with 
phononic gap. We consider a phononic crystal formed by periodic holes in Ge plate or a 
rigid cover / Ge layer / rigid substrate structure with gaps ~ a few GHz. The spin relaxation 
is suppressed dramatically, if the Zeeman frequency uz is in the phononic gap, but an effec- 
tive coupling between the spins of remote donors via exchange of virtual phonons remains 
essential. \iu>z approaches to a gap edge in these structures, a long-range (limited by detun- 
ing of uz) resonant exchange interaction takes place. We estimate that ratio of the exchange 
integral to the longitudinal relaxation rate exeeds 10^ and lateral scale of resonant exchange 
~0.1 mm. The exchange contribution can be verified under microwave pumping through os- 
cillations of spin echo signal or through the differential absorption measurements. Efficient 
manipulation of spins due to the Rabi oscillations opens a new way for quantum information 
applications. 


1 



Introduction 


Successful implementation of quantum information processing (QIP) requires not only invention 
of new quantum algorithms such as Shor algorithm, quantum error correction code or quantum 
adiabatic algorithm, but also further hardware development, i.e. realization of various qubit ar- 
chitectures - from trapped atoms to superconducting circuits.^ A significant advantage of solid 
state systems, based on different types of quantum dots^ or impurities in semiconductors,^"^ is a 
capability to fabricate, manipulate and read out qubits using semiconductor nanotechnology and 
conventional electronics. On the other hand, all reliable and efficient QIP schemes simultaneously 
require both long qubit decoherence times and controllable qubit manipulation, which poses a ma- 
jor challenge for practical implementation of these systems.® Indeed, shallow donor spin qubits in 
semiconductors have a number of advantages related to these requirements due to tunable spin- 
lattice interaction and a possibility to control spin states without a charge-induced noise. At the 
same time, broadly investigated silicon-based donor qubits with large spin dechorence times suffer 
from limitations in controlling and manipulating spins due to weak spin-orbit interaction.® Here 
we suggest a route for implementing new spin-qubit architectures based on donor spins embedded 
in specially crafted germanium structures (quasi-two-dimensional periodic phononic crystals or 
planar phonon waveguides) with large spin-orbit interaction of the Ge host and engineered phonon 
bangaps to simultaneously suppress spin decoherence and enable strong spin-spin coupling be- 
tween the qubits. 

The large spin-orbit coupling inherent to shallow donors in bulk Ge enhances our ability to 
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manipulate the spin qubits but could be detrimental for their coherence. In fact, the spin relaxation 
time of donors in Ge bulk is three to four orders of magnitude shorter than that in Si. To cope with 
this problem we will utilize Ge-based artificial periodic structures, known as phononic crystals 
(PC).^ Similarly to the photonic crystals, that were invented to control the light,^ the phononic 
crystals of different dimensions can be used to control various types of acoustic waves, e.g. to 
filter and focus sound ^ or even to create the earthquake proofing of buildings A proposed quasi- 
2D phononic crystal formed by a lattice of holes in a suspended Ge layer is shown in Fig 1-a. This 
structure is similar to recently manufactured Si-based PCs® with ~100 nm period and thickness 
defining the phonon gap within GHz frequency domain. The phonon dispersion curves, shown 
in Fig. lb, display pronounced gap in the frequency interval 13-^15 GHz. If the Zeeman energy 
of the donor spin fnoz is tuned inside the phonon gap the one-phonon spin-flip transitions will be 
forbidden due to energy conservation. As a result, the longitudinal relaxation rate, determined 
by very weak two-phonon processes, will be suppressed by five orders in magnitude compared to 
its bulk value. At the same time, the spin-lattice coupling will remain strong (3-^4 orders larger 
than in bulk silicon) and the spin-spin interaction via virtual one-phonon exchange processes will 
exceed by many orders. 

Similar behavior is possible if Ge layer is sandwiched between the rigid substrate and cover 
layers made from diamond, BN, or H-SiC (r/Ge/r-structure, see Fig. 2a). Penetration of vibrations 
from cover and substrate layers into Ge is weak because the reflection is effective due to the about 
10 times differences® between modules of elasticity in rigid materials and soft Ge. Thus, the spin- 
phonon interaction with bulk modes is ineffective and there is a gap for waveguide modes in Ge 
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Figure 1: Phononic crystal, (a) Geometry of PC formed by square lattice of holes, (b) Dispersion laws 
along rx, FM, and XM directions with gap between 6th and 7th modes, (c) Distributions of displacements 
over the unit cells of symmetric (left) and nonsymmetric (right) PCs for 7th mode at M-point (colored from 
blue, correspondent zero to red); sizes?dimentions are in nm. 


layer up to cut-off frequency in GHz range (Fig. 2b), so that is suppressed but the exchange 
interaction via vibrations localized in Ge remains essential if tiojz is in the phonon gap for waveg- 
uide modes. Once again, one obtains the donor spin system with an effective exchange interaction 
and a negligible vi. If ioz is fine tuned to lie within a gap edge in these structures, both strength 
and lateral scale of a resonant exchange interaction (REI) enhance. Thus, a spin system with es- 
sential?strong (and long-range for REI-regime) interaction and suppressed relaxation (transverse 
rate i >2 remains only) in PC or r/Ge/r-structure. 


Both verification of the spin-Hamiltonian parameters and manipulation of spin coherence 
are possible under resonant microwave (mw) pumping of frequency u u>z. Under continuous 
mw pump, the exchange-renormalized coz modifies a differential absorption shape in the linear 
response regime. An interplay between the Rabi oscillations frequency cor and the exchange con- 
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Wavemimber 

Figure 2: Sandwiched Ge structure, (a) r/Ge/r structure with Ge layer of width d\ incident {i-), 
reflected (r-), and transmitted (t-) waves and three waveguide modes (shear and coupled) with in- 
plane wave vector q are shown, (b) Dispersion laws (frequencies in 2^/qq/c? versus wavevector 
in 2/d) with gap up to cut-off frequency irct/d. 

tributions to Zeeman frequency (cji 2 — coz, see Fig. 3a) takes place for the nonlinear pumping case 
and it modifies a nonlinear differential absorption. If dephasing processes and long-range disorder 
are essential, a two-pulse spin echo measurements^^ enable to extract the exchange renormaliza- 
tion of u>z- Fig- 3b shows the sequence (vr/d — r — 37r/4 — r ^ echo signal) with different 
frequencies of free rotation during r-delay intervals. As a result, an echo amplitude oscillates with 
T, if C(Ji 2 7^ coz- Beside of this, a multi-pulse spin echo scheme can be applied for manipulation of 
averaged spin orientation. 


Summarizing the results obtained, we have demonstrated that a controllable manipulation of 
in quasi-free spin system in PC or r/Ge/r-structures is possible by the reasons: 
i/ strong suppression of relaxation-to-exchange ratio (in contrast to the bulk case^^) opens a way 
for fault-tolerant operations; 

ii/ very sharp (at deturnings Sz < 10“^ corrrespondent a weak variations of magnetic field) trans- 
formation from free spin system to large-scale REI-regime permits a remote control of qbits, with- 
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Figure 3: Spin dynamics under mw pumping, (a) Bloch sphere with Zeeman frequencies 
renormalized due to exchange, uji < uz < 002 , and Rabi oscillations frequency, ur, here 
Sj =const due to the spin conservation law. (b) Spin echo (7 t/ 4 — r — 37 t/ 4 — r) with 

r-dependent amplitude of echo signal due to difference of frequencies coi and U 2 . 

out any electric circuits when noise is suppressed; 

iii/ effective control of spin conversion between St^± and by microwave pulse takes place due 
to the interplay between exchange and Rabi oscillations; 

iv/ formation of macroscopic spin patterns, with lateral sizes up to mm, using micromagnets in 
order to control REI-regime, see and references therein; 

v/ manipulation of single spins (or a few-spin clusters with exchange interaction) in ultra-pure Ge, 
with inter-donor distances >3 fxm (concentrations < 10^^ cm“^), employing STM-technique, 
see^° and references therein. 


Results 

Exchange via phonon modes. The spin subsystem of donors in Ge is described by the Hamilto- 
nian Hg = (^/2) with the Zeeman frequency ojzHOZ if magnetic field is perpendicular 
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to the Ge layer and is the Pauli matrix for kth donor placed at r^. Within the harmonic approach, 
the quantized acoustic modes is described by Hph = and the spin-lattice interac- 

tion is introduced in analogy to the bulk case,^^ taking into account modifications of the vibration 
spectra in PC or r/Ge/r structure (see details in Supplementary Note 1 and Supplementary Figs. 
1-4): 

Hsph = 9uq^k • + H.C. (1) 

i/q 

Here 6+q and are the creation and annihilation boson operators for mode of frequency o;,^q 
with 2D wave vector q (which varies over the first Brillouin zone in the case of PC, |q| < 7r/a 
for square lattice of period a) and u labeled polarization and discrete quantum numbers. The 
components of vector are expressed through the strain tensor for donor at (averaging over 
the donor volume as provides the cut-off for wavelengths shorter the effective Borh radius^^’^^ , 
as, see Supplementary Note 2). The coupling constant p^q = KhuJz\Jb,/2puJi,cy is determined 
through the Ge density, p, the normalization volume, V, and the dimensionless factor K ~(1 or 
7.5) X 10^ for As or Sb donors. 

Within the second order accuracy with respect to Hspk> we transform the total Hamiltonian 

into^^ 

ffsphit) = E + , 

lyqka 

( 2 ) 

where the circular components of spin are dka{t) = (^ka ex.p{iau) zt) with a = ±1 and the plane- 
wave factors exp(iq • x^) are separated. After averaging of Heff over the phonon thermostat 


^e// — Hg+Hph-\- 


u 
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with the density matrix fjph (Tr^/j ... is the trace over boson variables) we obtain the effective spin 
Hamiltonian for the donor spins, H^fj = Tvph7]ph ((Heff ~ Hph^, which takes into account the 
renormalization of Zeeman splitting (polaron effect) and the exchange interaction via phonons, 

%s. 


The polaron contribution is added to Hg which contains now the renormalized Zeeman fre- 
quency (jJzk = coz + Awfc. The ratio Auk/coz is up to 2 x 10“^ for Sb donor and is weakly 
dependent on temperature and on donor position in unit cell of PC or across r/Ge/r-structure, see 
details in Supplementary Note 3. The exchange Hamiltonian is given by 

HsS — XI Jaa'^^ka^k'a', (3) 

aa'kk' 

{k^k') 

(kk') 

where the exchange matrix ^ has a non-zero circular components only. This matrix is given by 

2 / A^<i) fAci)* 

j{kk') _ %q j Srfe/,-a'Srfca _|_ ‘^^k-a^rya' ^ig.Axt.t./ 

\ ^I'ci + Oi'lOz - Ol'uJz 

and there is a strong large-scale dependence on the inter-donor distance Ax^^/ = — x^/, while 

dependencies of on are vanished under averaging over transverse coordinates z^' [here 
r = (x, z)] and over lateral donor positions in the unit cell of PC. 

Donor-based qubit. We consider the averaged exchange integral replacing the exponen- 

tial contribution (exp(iq • Ax^^/) by the zero order Bessel function Jo(gAx) where Ax is the 
mean inter-donor distance. After replacing — (^z) where c a character- 
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istic velocity and coefficients ~ 1, the ratio of the exchange integral to the Zeeman energy is 
estimated as 


J /\x ^ 



hcOz 4:7Tpd(P 


Qm 

dqqJo(qAx) 


K^hwz 




A 


Axi 


Trias. 


(5) 


Here we introduce the factor A^x = lo’^ dq . . . and Ma* ~ pdAx'^ estimated mass of 

Ge vibrated between spins. Here A/^x ^ 1 because the cut-off factor ^J^lAx/aB in /J- dq... 
and J^x decreases slower Ax~^. We perform numerical estimates of (5) for Sb donor assuming 
d czlOO nm. Using a;^ =15 GHz (if uz is approaching to the gap edges, an additional resonant 
enhancement of A/^x takes place) and taking c as an averaged sound velocity of Ge we obtain 
JAx/frcoz ~ 3t4ax X 10“^ i.e Jax/^ ranges up to 30-^90 kHz if Ax ~ d. 


This result should be modified under the REI-regime, when uz approaches edge of gap, cog, 
and Ax > d, so that the near-edge mode gives a main contribution to Jax- For these conditions, 
Jax/^z is transformed into 


Jax K'^huz T dqqJo{qAx) K'^huz 


hujz ^TTpd 2col J Sl + {q/qiy MaiduzY 


(2d'' 

Ko (SzqiAx ) , M, = p - 

V 7T , 


( 6 ) 


where the energy Md{duJzY in the denominator is estimated with the use |^ol ~ tt/o? and parabolic 
dispersion of mode, ujg[^ + {q/Qi)\ written through qi ~ ix/d. The argument of the modified 
Bessel function of second kind, Kq{ii), is determined by the relative deturning 5z = \J {oig — rcz')!oJz 
multiplied by Ax. The exponential suppression of Jax for Ax > dji^bz and slow (logarithmic) 
variations of Jax if c? < Ax < d/irSz take place due to the asymptotics Ko{y) cz ln(2/y) — C or 
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Y^7r/2y exp(— y) if y > 1 or |/ < 1; here C is the Euler’s constant. For the parameters used, Eq. (6) 
gives J/^x/'huJz ^ 8.8ii'o(. . .) x 10“® and Jax/^ is in agreement with estimates for (5), if /\x ~ d. 

In Fig. 4 we show the ratio Jax/^z for the sandwiched r/Ge/r-structure versus inter-donor 
distance, which demonstrates a sharp suppression of the RIE-regime if Sz increases, see inset. If 
7t5z > (disorder effects are discussed below), the REI-regime is realized over a wide 

interval of A.x, between 10 nm and 0.01^0.1 cm. Similar dependencies of Jax for PC with the 
dispersion relations of Fig. lb show the REI-regimes if uz approaches to the upper or lower edges 
of gap (M or T points of the Brillouin zone). The results for PC are more sensitive to lateral 
and transverse positions of donors in unit cell and, after averaging over donor distribution, J^x 
differs from the above estimates by factor 0.5-^3. In general, the the exchange matrix in (3) can be 
modified using selective doping (lateral or across Ge layer) or orientation of magnetic field with 
respect to OZ or to axis of PC. 

Further, we turn to the longitudinal relaxation rate z/i in the structures under consideration. 
In PC at low (<C fuvz) temperatures, ui is determined by the two-phonon spin-flip transitions while 
in the r/Ge/r-structures ui is caused by the spin-flip transitions via bulk phonons weakly penetrated 
into Ge layer, see Supplementary Note 4. We estimate the relative relaxation rate based on the 
golden rule with the forth-order (oc K"^) matrix elements of in PC and accounting bulk modes. 
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Figure 4: Exchange integral. Ja* versus Ax (in units 10~^hu)z and d, respectively) for Ge layer 
sandwiched between rigid slabs under relative detunings 10“^ (1), 5xl0“^(2), 0.05 (3), and 
0.5 (4). Asymptotics Jax oc Ax“^-® is shown by dotted curve and inset demonstrate In-dependency 
of Jax/^z on dz for irAx/d ~5 and 10 (solid and dashed curves, respectively). 


weakly propagated through r/Ge/r-structure. The results are: 


^1 

< 

UJZ 


{K^hujz/Mic^)\ 

K'^nuz/Mii^, 


Mi ^ pd{c/cuz)^ 
Mu = 2T^p{c/uzf 


(7) 


Here the characteristic sound velocity c is combined from q ^ in Ge and rigid materials and oc cr 
in r/Ge/r-structures. For the above parameters, we got vxjbJz < 10“^^ and z/f ^ exceeds seconds. 
In PC Pi oc uj\ nnd longitudinal relaxation time i/f ^ exceeds hundreds seconds. Thus, Jax/^^i 
exeeds 10^ and in photonic gap P\ is negligible. 


Bloch equation. Temporal evolution of weakly interacted spins randomly placed in PC or r/Ge/r- 
structure is described by the averaged spin vectors,^^ su = {l/2)Tvso-kPt, where pt is the multi- 
spin density matrix governed by the equation with the Hamiltonian 77e//- We restrict ourselves 
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by the second order accuracy on Hss (the mean field approximation) and factorize the two-spin 
correlation function as Tr sPt^k,ad'k',i 3 ~ ■Skt,aSk't,p, cf. Supplementary Note 5. As a result, the 
system of the nonlinear Bloch equations for Skt takes form: 

^ + 7 • Sfct = [{<^zk + AoJt) X Skt] - I • Sfc't) X Skt] ■ (8) 

™ ^ k'ik'j^k) 

Here A<x>t is the time-dependent Zeeman frequency under a microwave pumping (below Au;^ _L0Z), 
7 • Skt describes relaxation of kth spin, and the last term is written through the effective exchange 
frequency with matrix (4). Under pumping Aojzt switched on at f = 0 Eq. (5) should be solved 
with the initial conditions Skt=o = (0,0, s^o) where Sko — —1/2 if temperature <C hu>z and 
■Sko — ^ 0 for high temperatures (S> hujzk)- 

Instead of a microscopic set we introduce the averaged spin orientation 

Sxt = {J2k / {T.k ’ 

which is weakly dependent on transverse coordinate 2 ;. For the large (:^ d, a) scale inhomogeneity 
case Sxt is governed by the spin diffusion equation 

^ + 7 • Sxi = [Oxi X Sxi] - D [(AxSxi,x) X Sxt] , 

^xt = ^Zx — AxSxt,± + AuJt , (9) 

where uJzx takes into account the non-uniform Lamb renormalization of uz and the exchange 
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contribution is transformed into — WxSxt,±, where Sxt,± is the transverse part of spin orientation. 
The frequency uix is determined by the averaged exchange integral with the non-zero 

and equal xx- and r/r/-components. The diffusion coefficient, D, is estimated as ai/ex/2 where l^x 
estimates a scale of exchange interaction. The last contribution of (8) is negligible for the case of 
large-scale lex) inhomogeneities of ujzx and cDx- Because the ratio ^z^i/Jax is negligible, we 
replace 7 • Sxt by V2^xt,± with the transverse relaxation rate U2. 

Taking into account U 2 and neglecting diffusion if scale of disorder > I ex, one obtains the 
nonlinear (with respect to and Auzt) system for the transverse and longitudinal (e^^'x*,!!) parts 
of spin orientation 


+ ^2^ Sxt,± — ^xt,\\ [Oz X Sxt,J_] — A.U>tSxt,\\ , 

^5'xt,|| = (Au;* • Sxi,±) . (10) 

Here ^xt,\\ = <^zx + ^x>S'xt,|| includes the Lamb shift and the exchange (oc Sy^t,\\) renormalization. 
Within the collisionless regime, U 2 t <C 1, the spin conservation takes place y = 

with the xf-independent initial orientation Sq. If Acut — )■ 0 and 1^2 — ^ 0, Eq. (10) describes free 
rotation of Sxt,± around OZ with the frequency ujzx + cOxSq. Characterization of PC and r/Ge/r- 
structure (exchange, relaxation, and disorder parameters) and manipulation of spins are possible 
under resonant microwave pumping. 


13 



Microwave response. The absorbed power is given by Pxt = —hAujt • dSxt,±/dt, where (...) 
means averaging over period 27r/a;. Under weak circular pumping = copexp{iut), 

the solution of linear Eq. (10) gives the resonant peak = huuj,V 2 / [^2 + — <^x/2)^] where 

(iojx = cc — cc^x is the frequency detuning and Sq = —1/2 for zero-temperature limit. The resonant 
line has linewidth determined by 1^2 and by disordered contributions stem from Sojy^ and uix- If 
line is narrow enough, these contributions can be verified from the shape of dP/ dH averaged over 
disorder, similarly to the measurements of GeSi dots.^^ The derivative dP/dH increases under the 
REI-conditions due to additional dependency of oj on 5z- 

In the case of weak exchange, uj <C max(ccp, |(ic(j|, 1^1, 2), the linear (oc 5*0) system (10) de- 
scribes evolution of the resonant absorption Pxt and the spin orientation Sy^t,\\- Neglecting damping, 
at z/ 2 f <C 1 and z /2 <C ccp, |(5a;x|, and using the rotation wave approach, if \5uj\ ojz, one obtains 
oscillating responses 


Up 

Pxt Sq?IU sin CCxpf , pQ 00 s CCx7?f (11) 

UxR 

with the Rabi frequency cCxr = \/up + du^ and the tt/ 2 phase shift between P^t and Pxt,||- If the 
exchange interaction is essential (u ~ up, |5o;|), shape and strength of temporal Rabi oscillations 
are sensitive on u/up. Within the rotating- wave approximation, we plot these responses in Fig. 
5 for low temperatures, Pq = —1/2 at resonant condition (ia; = 0 (implicit solution for Sxt,\\ can 
be written through the elliptic integrals, see details in Supplementary Note 6). Shapes of temporal 
Rabi oscillations of Pj and St,\\ are modified with a suppression of their amplitudes and doubling of 
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periods around i^/ujp ~5.5 due to an effective interplay between exchange and Rabi oscillations. 
Detuning {5io ^ 0) results in a shift of similar behavior iouj/ojp ~8. The collisionless response 
is suppressed by decoherentization if V 2 t > so that an effective manipulation of Sj is possible 
with the use of pulse duration tp < ^2 • More sharp peculiarities take place for the differential 
nonlinear absorption dPt/dH. 



a>p t 


Figure 5: Interplay between exchange and Rabi oscillations. Absorbed power Pf (in units 
(jjpfiuj) and the longitudinal spin orientation 5'f,|| versus dimensionless time, copt, forCl/cop =2 (1), 
4 (2), 5 (3), 6 (4). 

The two-pulse spin echo scheme (7r/4 — r — 37r/4 — r echo) permits verification of 
exchange contribution under an essential decoherentization and long-range disorder. Here the 
delay times r ^ 2 and pulses of frequency u> and durations 2 are connected with pumping 

levels u)pi^2 by copiti = tt/ 4 and ujp2t2 = 37t/4. In the rotating-wave frame, free evolution 
of St,± = {St,x + iSt^y) exp{—icot) after first and second pulses is oc exp{Suji^2t) with different 
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frequencies Suji ^2 = — ^^<S'i^ 2 || and (. . .) stands for averaging over long-range disorder. If 

<^pi ,2 ^ ^2 and exchange is negligible according to Fig. 5, the spin orientations after first 

and second pulses are 5'i,|| = Sq/\/2 [c.f. Eq. (1 1)] and 5'2,|| = — 5'o(l -h cos 5u>it) /2. For the case 
of Gaussian disorder with the averaged variations of Zeeman frequency = 5u>*, spin echo 

signal is 


^t,± = -S'oexp 


[Su*{t — 2r)]^ 
2 


^ (0r) , 


( 12 ) 


where 5'o = (1 -I- \/2)/4 stands for the echo amplitude at t = 2r if a; — )• 0 and 4/(0) = 1. Function 
4/ {(pr) with (f)T. = SqU}t!\/2 describes modulation of caused by the exchange-induced differ- 
ence in 5uji and Suj 2 , which results in an interference oscillations of {(pr)- The disorder-induced 
exponent and the modulation 4/ are multiplied because of additive contributions of these factors to 
the frequency ^xt,\\ in Eq. (10). Shape of modulation of St= 2 r,± versus delay time r is determined 
by and Im4^ plotted in Fig. 6. Here and Im4^ are even and odd functions of pr and 
4>r > ^ corresponds to the spin inversion case, Sq > 0. Thus, verification of exchange contribution 
require variations of r in ~ 10 /iS scales. 


Discussion 

It is important to stress that these results are based on the estimates written through the ratios of 
K'^hiOz to the characteristic energies (factors M .c^ in denominators), which are evident from the 
dimensional requirements. The evaluation of Hgs is restricted by donor concentrations riD 10^® 
cm“^, when the exchange due to tunneling overlap of donors is negligible. For lower concentra- 
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Figure 6: Oscillations of spin echo amplitude. Real and imaginary parts of versus delay time 
T and 5*0 (here 0^- = j\f%. 


tions, the mean-field approximation for exchange in Eqs. (8-10) is valid if < radius of 

interaction. A region of intermediate concentrations, between the mean-field regime and a system 
of free donor spins (if < 10^^ cm“^) requires a special analysis. 


In order to implement the structures suggested and to measure the peculiarites found, ?? 
one have to meet several technological requirements for PC or r/Ge/r-structures: a) suppression of 
spin decoherentization rate V 2 in comparison to typical values in Ge-based materials^’^^ b) reduced 
stresses, dislocations, and interface disorder, c) homogeneity of donor distribution in PC or r/Ge/r- 
structures or possibility for selective doping, far from imperfections at boundaries or interfaces, 
and d) spatio-temporal stability of magnetic field and pumping characteristics, i.e. uz, cc, and ccp, 
as well as frequency of gap edge allowed realization of REI-regime. 


To conclude, we have demonstrated that a combination of phonon engineering provided gap 
in vibration spectra in PC or r/Ge/r-structure with unique spin and technological characteristics of 
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Ge opens a way for quantum information applications. We believe that our paper will stimulate 


effort for preparation of structures suggested and for verification of qbit parameters. 


Methods 

Theoretical background. Spin-lattice interaction in structures with phononic gap is considered 
by the classical proceedure^^’^^, based on the elastodynamics of continuous media.^^ The Bloch 
equation for averaged spin density is evaluated within the mean-field approach, which corresponds 
to the long-range exchange interaction. The nonlinear spin dynamics is examined with the simpli- 
fied description of decoherentization through V2 for the long-range inhomogeneities case. 

Numerical simulation. The eigenmodes and eigenstates for for the elastic vibration problem in 
PC or in r/Ge/r- structure are analyzed with the use of the finite elements method realized via Corn- 
sol Multiphysics software. The PC was designed as the structural mechanic modulus (unit cells) 
connected via the in-plane periodic (Bloch-Floquet) boundary conditions. The free boundary con- 
ditions are applied at all other surfaces. Standard numerical methods have been used for simulation 
of the exchange interaction and the nonlinear spin dynamics (the interplay of Rabi oscillations and 
exchange, or the spin echo modulation). More details of our calculations can be found in Supple- 
mentary Materials. 
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